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Abstract 

Deviations from kinetic equilibrium of massive particles caused by the uni- 
\q [ verse expansion are calculated analytically in the Boltzmann approximation. 

For the case of an energy independent amplitude of elastic scattering, an exact 
partial differential equation is derived instead of the usual integro-differential 
' one. A simple perturbative solution of the former is found. For the case of an 

energy-dependent amplitude the problem cannot be reduced to the differential 
equation but the solution of the original integro-differential equation can be 
found in terms of the Taylor expansion, which in the case of aconstant ampli- 
tude shows a perfect agreement with the perturbative solution of the differential 



equation. Corrections to the spectrum of (possibly) massive tau-neutrinos are 
calculated. The method may be of more general interest and can be applied to 
the calculation of spectrum distortion in other (not necessarily cosmological) 
nonequilibrium processes. 



1 Introduction 

Abundance of stable relics from the hot epoch of the universe evolution are usually 
calculated with the help of the integrated in momentum Boltzmann kinetic equa- 
tion, so that the complicated integro-differential equation is reduced to an ordinary 
differential equation for the number density of particles in question. To make this 
reduction one has to assume that these particles are in kinetic equilibrium with the 
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cosmic plasma so that their distribution is given by the expression 

f(p,t)=expMt)-E)/T(t)] (1) 

where E is the particle energy, fi(t) is an effective chemical potential (in total equilib- 
rium usually fi(t) = 0) and T(t) is the plasma temperature. Under this assumption 
the particle number density 

n ®=9.Jj0jifM (2) 
was shown to satisfy the equation |2| 

n + 3Hn = (av) (n 2 eq - n 2 ) (3) 

where H = a/ a is the Hubble parameter, (av) is the thermally averaged product 
of the annihilation cross-section and the particle velocity, and n 2 is the equilibrium 
number density obtained from the expression (^j) with // = 0. This equation was used 
in particular for the derivation of the cosmological bound on the mass of a possible 
heavy neutrino [|, f| and for the calculation of the mass density of heavy particles, e.g. 
stable supersymmetric relics, which may constitute invisible matter in the universe. 
A detailed derivation of eq. (|3|) valid for any form of the cross-section is given in 
ref.0 (see also the book 0). 

Usually kinetic equilibrium is maintained long after breaking of chemical equilib- 
rium, characterized by a nonzero //. Though elastic and annihilation cross-sections 
are typically of comparable magnitude, the rate of annihilation, (n/n) ann ~ avn, ex- 
ponentially goes down with decreasing temperature together with the number density 
of the nonrelativistic annihilating particles, n ~ exp(— m/T), while the probability of 
elastic scattering, (n/n) e i ~ avn , dies down much slower, because no ~ T 3 . Here no 
is the number density of massless (or light particles) which elastically scatter on the 
heavy ones and preserve kinetic equilibrium. 
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However for the particular case of tau-neutrino with the mass in MeV range, 
deviations from the kinetic equilibrium may be essential and can change the results 
of the calculations |7|, || based on equation (|3|). Correspondingly the nucleosynthesis 
bounds on a possible mass of nu-tau would be different. Due to a nonzero mass, v T 
would cool faster in the course of expansion, than massless particles. This in turn 
would result in less efficient annihilation of v T v T and in their larger frozen number 
density. Thus the nucleosynthesis bound on the mass of v T would be stronger. In 
ref.jSJ a simple attempt to take into account nonequilibrium corrections to the v T - 
distribution in energy was made. Namely it was assumed that the equilibrium form 
of nu-tau spectrum is maintained down to a certain temperature T e i, determined by 
the strength of the elastic scattering, and below T e \ the distribution function is given 



by the expression valid for noninteracting particles, / ~ expfy (m 2 /Tj + p 2 /T 2 )]. 
The effect was found to be relatively weak. 

An important influence on the primordial nucleosynthesis may produce nonequi- 
librium electronic neutrinos originating from the interaction with nonequilibrium tau- 
neutrinos. There are two possible effects. First, there should be an excessive cooling 
of massless neutrinos due to elastic scattering on colder-than-equilibrium i/ r 's. This 
makes the production of v T in the collision of the massless neutrinos (inverse anni- 
hilation) less efficient and correspondingly leads to a smaller number density of v T . 
It weakens the nucleosynthesis bound. Another effect is the distortion of v e spec- 
trum due to energetic electronic neutrinos coming from the process v T v T — > f e ^ e . 
Such nonequilibrium i/ e 's give rise to a larger neutron-to-proton ratio. This effect was 
found to be stronger than the first one and an overall stronger bound on m Vr was 
obtained 0. 

In ref . [|1(| it was assumed that massless neutrinos maintain their canonical distri- 



bution in momentum but, as the massive ones, acquire a nonzero chemical potentials 
(see eq. ([[])). This assumption permits to reduce the problem to the solution of a 



system of ordinary differential equations. It considerably simplifies the calculations. 
Such overall heating would result in a weaker nucleosynthesis bound but the neglected 
spectrum distortion have a stronger impact on the nucleosynthesis |1. 



It was claimed in ref.||ll|, on the basis of numerical solution of the complete 
system of the exact integro-differential Boltzmann equations, that deviations from 
kinetic equilibrium generated by massive v T are very strong and this gives rise to a 
much weaker bound on the nu-tau mass than it was previously stated]?], ||. In view of 
the arguments presented above the sign of the effect should be just opposite. Recently 
the authors found an error in their computations (Madsen, private communication) 
and the result seems to be close to those existing in the literature. 

Because of the complexity of numerical calculations it may be interesting to find 
an approximate analytical way to calculate spectral distortion for both massive and 
massless neutrinos. In this paper the spectrum of massive particles, which have elastic 
interaction with the equilibrium bath of massless ones, is found. In the following 
section the partial differential equation for the spectral density is derived from integro- 
differential Boltzmann equation for the case of the energy independent amplitude of 
elastic scattering. This equation permits a simple perturbative solution. For the case 
of the energy dependent amplitude such simplification was not found but one can 
solve the integro-differential equation expanding the distortion into the Taylor series. 
For the case of a constant amplitude of elastic scattering, comparison of this method 
with the perturbative solution of the differential equation, derived in Sec. 2, is made 
in Sec. 3. The agreement is found to be very good. More practically interesting case of 
weak interactions with the amplitude proportional to particle energies is considered 
in Sec. 4. The results are applied to the calculations of the corrections to the spectrum 
of massive tau-neutrinos. 
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2 Differential kinetic equation for a constant am- 
plitude of elastic scattering. 



Let us consider elastic scattering of massive particles with mass m, momentum p, and 
energy E on massless ones with momentum k and energy uj. The basic Boltzmann 
equation in the Robertson- Walker expanding universe and in the limit of Boltzmann 
statistics has the following form: 

\Up)h(k) - f m is')h(k')] (4) 

Massless particles are assumed to be in thermal equilibrium so that fo(u) = exp(— oj/T). 
In the case considered in this section, the amplitude of elastic scattering Aq does not 
depend on the momenta of the colliding particles. All the integration but one in this 
expression can be explicitly done (see Appendix). Introducing new variables x = m/T 
and y = p/T and the new unknown function C(x,y) = expire 2 + y 2 )f m (x, y) and 
assuming that the temperature drops with the expansion in accordance with the usual 
law, T = —HT, we can rewrite this equation in the form: 

-\A \ 2 )me- u / 2 r°° dy'y f e~ u/2 



oAtt^uxv Jo 



QAn^uxy Jo u' 
[C(x, y) - C(x, y')\ [ e H^'l/ 2 - e {y+y ' )/2 ] (5) 

where u = y/x 2 + y 2 and u' = \/x 2 + y' 2 . 

Successively multiplying eq.(^) by exp[(u + y)/2] and by exp(— y) and each time 
differentiating in y we arrive after some straightforward algebra to the following dif- 
ferential equation: 

JC" + 2J'C = - 64 ff X V 2 a \e~vdy \eW 2 uyd x (Ce^)] } (6) 
\A \ 2 m 11 v /J J 

where prime means differentiation with respect to y and J is given by the expression: 
J(x, V) = > 2 r dze-P U-*)- \e- yl2 r dze-* ( 1 - (7) 

£ Ju+y \ Z I A Ju—y \ Z I 



It can be checked that J satisfies the differential equation: 

J" - J I '4 = -ye- u/2 /u (8) 

We will solve eq.(||) perturbatively, looking for the solution in the form C(x,y) = 
Co(x) + Ci(x,y) where Co (a;) is the equilibrium solution; it can be found from the 
condition of particle number conservation: 



C (x) / dyy 2 e- u{x > y) ~ 1/x 3 
Jo 



(9) 



If the rate of the elastic scattering is large in comparison with the rate of the universe 
expansion, or in other words, the dimensionless coefficient in front of the r.h.s. of 
eq.(|6]) is sufficiently small we can neglect C\ in the r.h.s. and explicitly solve this 
equation in quadratures. Twice integrating by parts we get: 



C[(x,y) 



6A7i 3 C (x) Hx 2 [e 



-u(x,y)/2 



m 



J 



9 x Cq ( , . y 
u{x,y) + 



C n 



u(x,y) 



2 
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X 
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Co u(x, z) 



(10) 



The function C\ (x, y) is found by the integration of this expression in y and the 
integration constant (which is a function of x) should be chosen from the condition 
of the particle number conservation: 



Jdyy 2 e- u ^C 1 (x,y)=0 



(11) 



In a simple case of the universe dominated by relativistic particles the Hubble 
parameter is given by the expression 

/ 4-7r 3 g* m 2 



H 



45 mpix 2 



(12) 
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where mpi is the Planck mass and g* is the number of relativistic degrees of freedom 
contributing into the cosmological energy density. Thus Hx 2 is x-independent and it 
is convenient to consider the relative quantity r f = (C' l /Co)(m/Hx 2 )(\Ao\ 2 /647r 3 ). It 
is presented in fig.l for x = 10 as a function of y and in fig. 2 for y = 3 as a function 
of x. With rising x the correction C\(x, y) is getting bigger; asymptotically for large 
x, C\ ~ x. We can trust these results till C\/Cq <C 1 at least in the region of y where 
the supression due to the factor exp[— u(x, y)] is not too strong, i.e. for y 2 ~ 2x. 

3 Power series expansion for the case of a constant 
amplitude. 

Here we will present an alternative approach to the solution of the original kinetic 
equation (||). This equation is not as easy to solve perturbatively (for small deviations 
from kinetic equilibrium) as the equivalent differential equation @. However the 
method considered here can be useful for the solution of the kinetic equation for the 
case of weak interactions when the amplitude strongly depends on the energies of 
the colliding particles and the reduction to the more simple form (Q) is not found. 
Below we will find the solution of eq.((^) expanding C(x,y) in powers of y (in fact of 
y 2 ). Formally the expansion is valid even if the elastic scattering is relatively weak 
and the deviations from the kinetic equilibrium are large, but the convergence of the 
expansion is getting worse for larger x, when the probability of elastic scattering goes 
down. 

After a simple rearrangement we can rewrite eq.([|) in the form: 




(13) 
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We will look for the solution of this equation in the form of the Taylor expansion 

d(x, y) = C w (x) + C 12 (x)y 2 /2 + C u (x)y A /A + ... (14) 



The coefficient C\q is not determined by the equation but should be found from the 
condition of conservation of the particle number as is discussed above. One can easily 
check that the odd powers of y in this expansion are absent. This was evident if fact 
from the very beginning because the expansion in terms of the three momentum |p| 
normally goes in terms of p 2 and there is no reason to expect a singularity ^[p* 1 in 
the considered problem. 

We will retain in C only terms up to the fourth power in y, so the last integral 
in equation ([14]) which contributes in the order y 5 can be neglected. After a simple 
algebra the following system of equations, determining C\ 2 and Cu, can be obtained: 

\c A + \ Cih = X {?f-,) (15) 



where 



Ux) = I ' ^ e -K-^')/2 (17) 
u 



For example for x = 10 we find C 2 = — 1 and C4 = 0.007. It very well agrees 
with the results of the previous section up to y ~ 10. At these and larger y's the 
suppression of the spectrum due to the factor exp(— u) is already very strong and 
this part of the spectrum is not physically significant. It is noteworthy that there is 
no formal small parameter in the equation (|15|) and all C n would in principle equally 
contribute. However higher order terms are small numerically, possibly because of 
nonrelativistic nature of the problem, and thus the neglect of those terms is justified. 
For example if we neglect the term C4?/ 4 (which we did not) the resulting value of C 2 
would be only slightly different form the one found together with C 4 . 



4 Weak interaction scattering. 

Now let us consider a more realistic case of elastic scattering of Dirac massive tau- 
neutrino on massless z/ e , u^, and e ± , which are in thermal equilibrium. The amplitude 
of elastic scattering summed over all light lepton species is 

\A W \ 2 = K[(pk) 2 + (pk') 2 ]/m 4 (18) 

where (for sin 6^ = 0.23) 

K = 32G|m 4 (3 - 4 sin 2 6 W + 8 sin 4 6 W ) « 1.09 x 10" 20 '{m/MeV) 4 (19) 

We have to substitute the amplitude A w instead of A into eq.(£|). The final result 
of the integration is grossly simplified in nonrelativistic limit when one neglects the 
terms of the order of p 2 /m 2 in the amplitude. The details of the integration can 
be found in the Appendix. Performing it we find, instead of eq.(H), the following 
equation: 



HdWOO = - 12 £TT r d -^*-^\c{x,y)-CM)] 

b / i'K 6 yx t> Jo u 

2 sinh(y/2)(l + y'/3 + (y 2 + y l2 )/2A) - 2 2 /cosh(y/2)(l/3 + y' /12)] (20) 

As in the previous sections we make the perturbative expansion C(x,y) = Cq(x) + 
Ci(x, y) and expand C\ in powers of y. In eq. (p0|) we have kept only the terms which 
are essential for calculations of C\ up to y 4 . The exact r.h.s. of this equation is 
presented in the Appendix. 

We will slightly change notations here and expand the relative deviation from the 
equilibrium distribution: 

r w (x, y) = C a (z, y)/C (x) = r 2 (x)y 2 /2 + r A (x)y 4 /A (21) 

For the coefficients r n the following equations are valid: 

i, f3 + i riFs = s m-i) (22) 



Here S = l67r 3 (Hx 2 )x A / Km A and 

F n = r dyy n - + % + e ( "^ +y)/2 (24) 
Jo u \ 2 8 J 

L n = [°° dyy n - e («-*+«)/2 (25) 

Jo u v 48 192/ v ; 

For example for m = 20 MeV and x = 10 we get r 2 = —4.5 x 10~ 2 and r 4 = 
6.66 xlO -4 . Behavior of r w (x, y) for x = 10 in units S/x A ~ 18.6Mev 3 /m 3 is presented 
in fig.3. The coefficients r 2 (x) and r±(x) in the same units are presented in figs. 4 and 
fig. 5. The corrections quickly rise with x (remember the factor x 4 / S) but they remain 
reasonably small till the freezing of the annihilation of tau neutrinos. 

Comparing the annihilation rate for v T in kinetic equilibrium with the rate of 
annihilation of v T with the spectrum calculated here, we conclude that the frozen 
number density of nu-tau can be bigger by about 10% than that found in the standard 
case. However, calculating the spectrum, we have neglected the annihilation of v T v T 
and this process may have an essential influence on the spectrum and on the final 
result on the frozen abundance. We will take this into account in the subsequent 
work. 
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A Appendix. 

Integration in the r.h.s. of eq.([|) can be done in the following direct way. First 
using 5 3 -function we integrate over d 3 k so that k = k' — p + jf . Correspondingly 
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uj 2 = k 2 = uj' 2 — 2u'Rcos8k + R 2 where R = \p — j7\ and 9 is the angle between R 
and k . 

Since dcosO, which comes from d 3 k', can be whiten as 

dcos9k = —duuj/uj'R (Al) 

we can easily integrate over duj using 5(u — uj' + E — E'). The limits of the integration 
over uj' are determined by the condition | cos^l < 1 so we get 

I fk^ {p + k _ p ,_ k , ) = *L r dJ{ ^ hk {A2) 

J UJ UJ -ft Jul miri 

where uj m in = (R+E — E')/2 and the expression in the angle brackets means averaging 
over the polar angle of the vector k'. For a constant amplitude there is no dependence 
on this angle. 

In the remaining integration over d 3 p' we change the integration over the azimuthal 
angle to the integration over R in accordance with the relation dcos8 p = —RdR/pp'. 
Correspondingly we get 

where, as above the averaging over p'-polar angle is done. The limits of integration 
in R are R min = \p-p'\ and R max =p + p'. 

If the amplitude of the scattering is constant and if the light particle distribution 
functions are known, f = exp(—E/T), we can easily make all the integration but one 
(in dp') and obtain eq.(|5|). For the case of weak interactions the amplitude fll8|) is a 
function of momenta and the integration is more complicated. We will make it in the 
nonrelativistic limit, neglecting terms of the order of v 2 = p 2 /E 2 in the amplitude. In 
this limit \A W \ 2 = K(E/m) 2 ((uj 2 + uj' 2 )/m 2 . The absence of the angular dependence 
in the amplitude simplifies the calculations very much. Performing the integration 
in the way described above, we obtain for the r.h.s. of eq.([|) (with the substitution 
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Aq — > A w ) the following expression: 

128lT- i pm 4: JO E' J\p-p'\ JuJmin 

(26) 

where us = u' + E — E' . Making a trivial integration over u and neglecting terms 
proportional to (E — E') 2 ~ m -4 and also neglecting the integral similar to the last 
one in eq. ([13|) we arrive to eq . (|20|) . 
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Figure captions. 

Fig.l. Relative deviation from the thermal equilibrium of the normalized y- 
derivative of Ci(x,y) for the case of a constant amplitude of elastic scattering for 
fixed x — 10 and running y, r'(10, y). See the definition of r' at the end of Sec. 2. 

Fig. 2. The same as in Fig.l but for fixed y = 3 and running x, r'(x, 3). 

Fig. 3. The correction to the spetral density of nu-tau for x = 10, r w (10,y) in 
units S/x 4 (see the end of Sec. 4). 

Fig. 4. Taylor expansion coefficient r^x) in units S/x 2 . 

Fig. 5. Taylor expansion coefficient r 4 (x) in units S/x 2 . 
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